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Abstract: We consider the capture of dark matter in the Sun by inelastic scattering against 
nuclei as in the inelastic dark matter scenario. We show that, assuming a WIMP-nucleon 
cross-section of a n = 10 -40 cm 2 the resulting capture rate and density are sufficiently high 
so that current bounds on the muon neutrino flux from the Sun rule out any appreciable 
annihilation branching ratio of WIMPs into T^ -1- !^", Z°Z°, t + t~, ti and neutrinos. Slightly 
weaker bounds are also available for annihilations into bb and cc. Annihilations into lighter 
particles that may produce neutrinos, such as pions and kaons are unconstrained since 

those stop in the Sun before decaying. Interestingly enough, this is consistent with some 
recent proposals motivated by the PAMELA results for the annihilation of WIMPs into light 
bosons which subsequently decay predominantly into light leptons and pions. 



1. Introduction 



In the past twenty or so years, much effort has been devoted to probing galactic dark matter, 
(DM) beyond its gravitational effects, in an attempt to uncover what consitutes DM and 
how it interacts. Maybe the most promising candidate is a weakly interacting massive par- 
ticle (WIMP). Direct-detection searches look for WIMP-nucleus scattering in underground 
detectors and carry with them the hope of having laboratory control over DM. So far, most 
experiments have reported null results and placed very strong bounds on the elastic cross- 
section for WIMP scattering against nuclear matter. Notwithstanding these results, the 
DAMA collaboration announced an 8.3 a discovery in the annual modulations of nuclear 
recoil rate [1]. Unfortunately, if attributed to WIMP elastic scattering, such event rates are 
excluded by more than two orders of magnitude by the other experiments. An intriguing res- 
olution of this controversy is the proposal of "inelastic dark matter" (iDM) laid forth in Ref. 
[2]. If DM scattering off of nuclear matter is inelastic, with a transition to an excited state 
of DM roughly 100 keV above the ground state, then the incosistency between the different 
experimental results can be settled [3]. 

Alongside the direct detection efforts, there are also indirect searches looking for energetic 
neutrinos from WIMP annihilations in the Sun or the Earth or other energetic fluxes (photons, 
positrons, and antiprotons) from annihilations in the intragalactic medium. Very early on, in 
a seminal paper, Press and Spergel [4] considered the possibility of WIMPs' capture by the 
Sun and their computation was later refined and corrected by Gould [5]. In their work, and all 
subsequent work on the subject, capture was assumed to proceed through elastic collisions of 
WIMPs with matter. In this paper we consider the iDM scenario where capture proceeds via 
inelastic scattering. We compute the associated capture rate and discuss detection prospects 
in neutrino telescopes. 

It is important to realize that the computation of the neutrino flux, starting from WIMPs' 
capture and ending with their annihilations, involves several steps which are model dependent. 
The problem nicely divides into a question of capture and a question of annihilation. These 
are logically disjoint and call for a modular approach which we carefully review in section 3. 
Throughout, we strive to maintain a model independent approach and clearly state where 
and which assumptions are made. 

In section 2 we make a short detour and consider the relation of the results presented in 
this paper to some of the recent developments associated with DM observations and model 
building. Section 3 is devoted to the formulation of the problem and the modular approach 
we adopt in resolving it. In section 4 we present and discuss the results and section 5 contains 
our conclusions. 

Finally, we note that in their original work on iDM, Smith and Weiner already considered 
the possibility of capture of WIMPs in the Sun. Our quantitative conclusions differ from theirs 
on three accounts. First, we include form-factor effects which are extremely important for 
iron. In the elastic case, these effects are so significant as to dethrone iron from its place as 
the prime capturer of WIMPs. Second, we take into account the case where elastic scattering 
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is altogether absent and show that the resulting density is nevertheless sufficiently high to 
saturate the annihilation rate in the Sun. This ensures maximal signal and closes a loophole 
whereby iDM models with little or no elastic component can escape the bounds presented 
below. Third, the oscillations of the different neutrino flavors into each other are now well 
established. This closes another loophole whereby sneutrino type WIMPs annihilate into 
electron-neutrino only. Previously, the possibility existed for these neutrinos to keep their 
flavor and therefore result in no upward going muons in Earth-bound detectors. 

2. Recent Developments 

While this work deals specifically with the problem of inelastic capture of WIMPs by the Sun, 
we momentarily diverge off this route to discuss its relation to some of the recent developments 
related to DM. This discussion is useful on its own right because it serves to illustrate the 
different ingredients that enter the computation and expounded upon in later sections. 

Motivated by the positron excess seen by PAMELA [6], a new class of DM models was 
suggested in Ref. [7], employing a light mediator that both explains the large cross-section 
and the absence of any excess in the anti-proton flux [8]. A very reasonable candidate for 
this mediator is the gauge-boson of an additional (non) abelian gauge group under which 
DM is charged. These models have the added benefit that the small splitting (~ 100 keV) 
and inelastic coupling associated with iDM are naturally generated 1 . These models form 
an organic whole, capable of explaining many of the astrophysical anomalies reported in 
the recent years, and contain several distinct phenomenological signatures important for this 
work, 

1 . Sommerfeld enhancement of the annihilations of slow WIMPs via a relatively light boson 
helps to reconcile the large rate needed to explain the PAMELA positron excess with 
the cross-section deduced from the thermal relic abundance of DM. 

2. WIMPs annihilate predominantly into a pair of these bosons which in turn are weakly 
mixed with the SM. Their decays into /U + /t/~ or e + e~ are responsible for the excesses seen 
in PAMELA, and possibly ATIC, whereas decays into anti-protons are kinematically 
suppressed [12]. 

3. The features seen in PAMELA (and possibly ATIC) require a fairly heavy WIMP, 
anywhere between 100 GeV — TeV. 

4. DM excited states are naturally present in the spectrum with 5 = m x '—m x ~ 0(100 keV— 
MeV) which can realize the iDM scenario. 

Nevertheless each of the above elements can independently arise in specific constructions. 
Hence, from a phenomenological point of view we should consider how each of these elements 

x The idea was first suggested in Ref. [7] for the non-abelian models and expanded upon in Ref. [9]. The 
existence of this splitting in abelian models was demonstrated in Refs. [10, 11]. 
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seperately impacts the prospects of indirect detection of WIMPs capture and annihilation in 
the Sun in neutrino telescopes experiments. 

Indeed, such a signal may be more robust than any of the other indirect probes. Anomalies 
in electromagnetic signals (PAMELA, ATIC) may be due to special astronomical mechanisms 
rather than WIMP annihilations. On the other hand a putative signal of energetic muon- 
neutrinos coming from the Sun during the (southern) winter or straight up from Earth's 
center may have no other reasonable explanation other than WIMP annihilations in the Sun 
or Earth 2 . 

The first ingredient involving the Sommerfeld enhanced annihilation cross-section is in 
fact a boon since it helps to gaurantee equilibrium as discussed in section 4. In particular it 
can greatly enhance the rate of neutrinos from the Earth as discussed in Ref. [14]. 

The second ingredient concerning the annihilation channels has the most dramatic effect 
as far as neutrino telescopes are concerned. If the WIMPs dominantly annihilate into light 
bosons which consequently decay into muons and electrons, then no energetic neutrinos will 
be observed. Electrons do not yield neutrinos, and the neutrinos from the muons' decay will 
be too soft since the muons first stop in the Sun and only then decay, yielding a Michel 
spectrum ~ 30 MeV neutrinos [15, 16]. So this class of models, in their purest form, lead 
to no observable consequences for neutrino telescopes 3 . In what follows, we will instead 
remain agnostic about the annihilation channels open for WIMPs and quote results based on 
branching fractions. 

The third ingredient pertaining to the mass of the WIMP is relevant for the capture 
rate (and only mildly to the equilibrium condition). The capture rate is reduced by roughly 
~ m~ 2 . One power of the mass is coming from the obvious inverse dependence of the WIMPs 
number density on their mass, once the local energy density is fixed, p x = 0.3 GeV cm~ 3 . 
A second power of the WIMPs mass is present because it becomes increasingly difficult to 
transfer enough momentum and gravitationally capture heavier WIMPs. The heaviest target 
in the Sun is iron so the capture rate for WIMPs with m x S> 52 GeV is strongly suppressed. 
The resonant enhancement in the Earth discussed by Gould [5] is altogether absent. 

This reduction in capture rate is, however, offset by the increased neutrino energy. First, 
the conversion probability of neutrinos into muons in the rock or ice below the detector grows 
as oiv^ + N — > \i + X) ~ E v . Second, the range of muons grows with their energy as well 
~ E^. This offset does not occur for annihilations of WIMPs in the Earth. The density of 
accumulated WIMPs in the Earth does not build to be "optically thick" to neutrinos and 

2 To put the issue of particle versus astrophysical explanation in perspective, recall the hotly debated 
" Solar Neutrino Anomaly" . Some of its experimental evidences were suspected and astrophysical explanations 
suggested. The late John Bahcall correctly argued that while some experiments suggesting neutrino oscilations 
for certain neutrino energies may have astrophysical explanations the sum total of all evidences cannot [13]. The 
present "Unified DM models" motivated by several anomalies in the electromagnetic spectrum from microwaves 
to multi-GeV Gammas may survive in some form or another. Note that also the solar case involved a more 
robust neutrino signal. 

3 This is also true for models where WIMPs dominantly annihilate into photons, gluons, and kaons, all of 
which result in little if any energetic neutrinos. 
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hence the signal is proportional to the the square of the capture rate and massive WIMPS 
yield a weaker signal. These considerations omit two other factors: i) The Sun becomes 
optically thick for a neutrino of energy E v > 1/2 TeV produced at the solar core, reducing 
the solar signal; ii) More energetic muons point more precisely in the direction of the parent 
neutrinos. The 10 -3 radian precision required to resolve the solar core is unavilable. However 
better pointing improves the Earth signal since the annihilations occur at the center and this 
can be well tested by the long vertical strings of photomultiplyers in IceCube. 

The fourth ingredient we comment on is the most relevant for the analysis in this paper, 
namely the inelastic transition. As clearly explained in [2], an inelastic transition of ~ 100 keV 
results in a preference towards heavier targets and is essentially how the absence of nuclear 
recoils in CDMS (germanium) can be reconciled with the DAMA (iodine) results. Since iron 
is lighter than germanium the rate for capture in the Earth is tiny at best (nickel is actually a 
little heavier than iron, but its abundance is two orders of magnitude smaller). The situation 
would have been just as bad in the Sun if it was not for the kinetic energy gained by the 
WIMP as it falls in the Sun's gravitational well. The escape velocity in the Sun ranges in 
v esc ~ 600 — 1300 km/s and provides sufficient energy to overcome the excitation barrier of 
~ 100 keV. This entire paper rests on this simple observation, one that was made already in 
the original work of Smith and Weiner [2]. 

Finally, we close this section with a brief survey of existing iDM models and their rele- 
vancy to this paper. Rather than concentrating on the DM identity, a more useful categoriza- 
tion is obtained by focusing on the mediator of the inelastic transition. One possible mediator 
is the SM Z° as happens for example in the case of sneutrino DM [2] (see [17, 18] for many 
other interesting possibilities). If kinematically allowed, DM can then annihilate to a pair of 
Z°'s, which as we show below is strongly constrained. Similar comments would apply for a 
heavier Z', however, in this case the annihilation may be kinematically suppressed. Mediation 
through a scalar (higgs or other) is somewhat less constrained since DM annihilation into this 
scalar may produce only bottom quarks as final products. As already mentioned above, the 
possibility of a new light mediator can avoid the bounds entirely. DM annihilation in the Sun 
into any light entity (< GeV) would result in no observable neutrinos. Other variants are 
certainly possible, and we hope that the modular approach of this paper will allow for direct 
comparison with any other future models of iDM. 

3. Methodology 

3.1 Approach to Equilibrium 

The time evolution of the WIMP population in the Sun is controled by (we follow the notation 
ofRef. [19]), 

N = C-C A N 2 (3.1) 

where C is the capture rate of WIMPs in the Sun and Ca is related to the annihilation rate 
as Fa = \CaN 2 . The evaluation of the capture rate, C, will be presented below. In this 
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subsection we concentrate on the estimation of Ca and the modification necessary for inelastic 
scattering. Its significance manifests itself by solving Eq. (3.1) for the annihilation rate at a 
given time, 

T A = ictanh 2 (t/r cq ) (3.2) 

where r cq = 1/ ' \JCCa is the time scale required to reach equilibrium between capture and 
annihilation. When r cq becomes comparable or larger than the age of the solar system t Q ~ 
1.5 x 10 17 s, the system has not yet reached equilibrium and the annihilation rate is strongly 
suppressed, 

T A ^ l -C{t/r c ,f (3.3) 

On the other hand, when t Q » r eq , the annihilation rate is saturated at Ta = \C. 

In the usual WIMP scenario with elastic scattering against nuclear matter, Ca is essen- 
tially determined by the temperature of the Sun and the annihilation cross-section. This 
is based on the assumption that after being gravitationally captured, WIMPs will continue 
to lose energy with every subsequent collision and reach thermal equilibrium in the Sun's 
core. Even with elastic cross-sections as low as <r n — 10 cm the number of collisions they 
undergo over the solar lifetime is ~ t Q (a n v)n Q = 3 x 10 7 , which is sufficient to thermalize. 
In that case, their number density is given by [20], 

n(r) = n e- m ^' T = n^ 1 ''^ (3.4) 

where no is their number density in the center of the Sun, <j){r) is the gravitational potential 
with respect to the core, and T is their temperature. The second equality holds if we assume 
a constant core density p and define the thermal radius, 

/ 3T \ 1/2 ( T \ 1/2 /100 GeV\ 1/2 , , 

where Rq is the Sun's radius and T = 1.2 keV is the Sun's core temperature. The annihilation 
rate per WIMP squared is given by, 

_ fd 3 r n(r) 2 {<r A v) {cr A v) . . 



(/Anfr)) 2 (2vr) 3 / 2 



r th 



where a a is the annihilation cross-section, and the last equality follows from Eq. (3.4). The 
ratio of solar age to equilibrium time is then given by, 

to = 1Q 3 ( C \ V 2 / (oav) \ 1/2 / 0.01 xR Q \ 3/2 

r cq Vl0 25 sec-V \3 x 10" 26 cm 3 sec" V V r th J { ' 

which implies that equilibrium has been reached long ago and we can expect the full signal 
from the Sun. 

The situation is more subtle in the iDM scenario. First, let us consider the case where 
a small elastic component exists in the scattering of WIMPs against nuclear targets. In this 
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case, while capture proceeds through an inelastic transition (with a large cross-section), the 
subsequent thermalization of the WIMP against the nucleons in the Sun can be due to its 
suppressed (a n < 1CT 43 cm 2 ) elastic scattering. In specific models of iDM, the elastic coupling 
is suppressed with respect to the inelastic one by ~ 10 -6 4 . That would result in an elastic 
cross-section (per nucleon) of a n ~ 10~ 52 cm 2 which is too small to bring the WIMPs into 
thermal equilibrium with the rest of the matter in the Sun. However, if we take the elastic 
cross-section to be a < 10~ 43 cm 2 (current bound from direct-detection), but a n > 10~ 47 cm 2 
the WIMP can undergo enough collisions to thermalize. In that case, the above estimates for 
the annihilation rate Ca still apply. 

As was recently clarified and emphasized in Ref. [21], even in the absence of any elastic 
coupling, the second term in the Born series yields an elastic transition. When the force 
mediator is very light < 50 MeV, the resulting elastic cross-section per nucleon is sufficiently 
large (~ 1CT 42 cm 2 ) to allow for thermalization. However, the cross-section drops rapidly as 
the mass of the mediator increases and becomes inefficient for a mass of ~ GeV. 

In order to close any possible loophole in the argument we computed the final WIMPs' 
density in the pure inelastic case where no elastic scattering is allowed. In order to facilitate 
the computation, we approximated the Sun's gravitational potential with an analytic form 
which allows for an exact solution of the oribts 5 , as explicated in a beautiful paper by Henon 
[22, 23]. We describe the details of the computation in appendix B. In Fig. 1 we depict 
the resulting density for a particular choice of parameters. The WIMPs undergo only few 
collisions before their kinetic enery anywhere in the orbit drops below the inelastic threshold. 
Therefore, it is difficult for them to shed off all their angular momentum which is the reason 
why the density vanishes as we approach the Sun's center. On the other hand, since most of 
the Sun's mass is concentrated in its core, most collisions actually happen in the inner radius 
(< 0.2-Rq) and the resulting orbits are to a large extent contained inside the Sun as shown 
in Fig. 7. As the WIMP gets heavier, it become harder for it to shed off kinetic energy in 
every collision, but since its initial kinetic energy is larger, it also undergoes more collisions. 
Due to the different approximations used in arriving at these results, they are probably only 
good to 20% accuracy. However, considering that tanh(t Q /r oq ) pa 1 for to/r oq S> 1, these 
uncertainties do not propagate into the annihilation rate and muon yield discussed below. 
We conclude that even in the worst case scenario where no elastic scattering is allowed, the 
annihilation rate still reaches equiliberium and is saturated by half the capture rate 6 . 



4 This is usually achieved by starting with a Dirac fermion (or a complex scalar) of mass ra x and adding a 
small Majorana mass m. That splits the fermion into its Majorana components and generates a mass splitting 
of m x ' — m x = 5 = m 2 /m x « 100 keV. Any vector current coupeled to the fermion will result in a dominantly 
inelastic coupling between % an d x' ■ A small elastic coupling of size 8/m x « 10~ 6 — 10~ 7 is also generated. 

5 We would like to extend our warm gratitude to D. Lynden-Bell for bringing this solution to our attention. 

6 We have assumed here that after its first excitation collision, the WIMP quickly deexcite to the ground 
state so any subsequent collision is also endothermic. If this is not the case, and the deexcitation time is very 
long (as recently discussed in Ref. [24, 21]), the situation is even more favorable since the WIMP can now 
scatter against all the nuclei in the Sun, including hydorgen and helium. 
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Figure 1: On the left pane we show the WIMP's density (normalized to unity) against the distance 
from the Sun's center for m x = 100 GeV and S = 100 keV (solid-black), S — 150 keV (dashed-blue), 
and S — 200 keV (dotdashed-red) . On the right pane we depict the ratio t Q /r cq as a function of the 
WIMP's mass for 5 = 100 keV (black-solid) and 5 = 150 keV for (<r A v) = 3 x 10~ 26 cm 3 sec" 1 and 
C = 10 25 sec" 1 as in Eq. (3.7). 



3.2 Capture Rate 

If equilibrium is reached then WIMPs' annihilation rate is saturated at half the capture rate 
Ta = \C The formulae and procedure related to the computation of the inelastic capture 
rate, C, are presented in appendix A. In here we discuss some of the qualitative points 
that arise and clarify the modifications involved. For concretness, we fix the inelastic energy 
threshold to be 5 = 100 keV whenever numerical values are used. 

As the WIMPs fall into the Sun, they get accelerated in the gravitational well until they 
meet a nucleus against which they collide. Their velocity at that point is given by, 

w {rf = u 2 + v{r) 2 (3.8) 

where u is their velocity at inifinity, determined by the velocity distribution in the halo and 
the Sun's motion through the galaxy, and v (r) is the escape velocity from the Sun at radius r. 
For a WIMP to undergo an inelastic transition against another nucleus of mass m N , the total 
kinetic energy in the center of mass frame must be greater than the inelasticity, 

-fj,w(r) >6 ii = * 3.9 

2 m N m x 

This condition is more easily satisfied for heavier nuclei as we show in Fig. 2 with a plot of 
K.E. for the relevant elements in the Sun as a function of the mass enclosed in a given shell 
(plotting against the radius is more intuitive, however, it can be misleading since the mass 
density is not uniform). We took the velocity at infinity to be the most probable speed in a 
Maxwell-Boltzmann distribution of velocities with a Sun rotation velocity of 220 km/s 7 . We 

7 This fixed choice of velocity was made only with regard to Fig. 2. In the rest of the paper we used the 
full Maxwell-Boltzmann distribution for the velocity at infinity as explained in the appendix. 
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see that iron is an effective scatterer throughout the Sun, but oxygen is only useful in the 
inner 50% of the Sun's mass and helium is altogether useless as the kinetic energy is never 
sufficient to overcome the inelastic barrier. Indeed, as we shall see in section 4, scattering 
off of iron yields the largest capture rate. This is in contradistinction with the elastic case, 
where iron places only third behind oxygen and helium because of the form-factor suppression 
discussed below [5]. To obtain the results presented in section 4 we included all the different 
elements persent in the Sun with their proper densities [25] . The heavier elements were taken 
from [26] and their radial profile was assumed to follow the mass profile of the sun. With 
the exception of hydrogen and helium, all the different elements included allow for scattering 
somewhere in the Sun and contribute to the capture rate. 




m (units of M SUTI ) m (units of M sun ) 



Figure 2: The kinetic energy in the WIMP-nucleus center of mass frame as a function of the mass 
m (in units of Mq) enclosed in a given shell. The velocity at infinity was taken to be u = 220 km/s 
as explained in the text. The curves correspond to scattering against iron, oxygen, and carbon from 
top to bottom. The WIMP mass was fixed at 100 GeV (500 GeV) on the left (right) pane. Carbon is 
the lightest element present in the Sun against which the WIMP can scatter inelastically. 

The inelastic barrier in fact enhances the capture rate in those cases where a scattering 
is kinematically accessible. In order to be gravitationally captured the WIMP must lose a 
certain amount of kinetic energy. In the elastic case, the entire energy must be transferred 
to the nucleus against which the WIMP scatters. In the case of iDM, however, a significant 
fraction of this energy is lost to the excitation. When the excited state relaxes back to the 
ground state, it can only do so through the emission of very light states with no significant 
recoil. The net result is that the inelastic transition allows the WIMP to lose energy efficiently 
and hence to further enhance the capture rate. 

The transfer of a large fraction of the kinetic energy to an excited state instead of to the 
nucleus further contributes to the rate by softening the momentum transfer and taming the 
reduction usually associated with the nuclear form-factor. This is particularly important in 
the case of iron where the form-factor suppression in the elastic case is very significant and 
can result in an order of magnitude reduction in rate for m x ~ 100 GeV. More specifically 
consider a form-factor exp(— Q/Qq) with Qq = 82 keV for iron (see appendix A for details) 
and where Q is the energy transfer. If Q is now smaller by about 5 ~ 100 keV this results in 
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a rate increase by a factor of exp(<5/Qo) = 3.4. 

The above considerations lead us to an important conclusion. While one might have 
naively expected that inelasticity will result in a substantial reduction in the capture rate 
this is in fact not the case! Iron (which in the elastic case would have dominated capture 
if not for the form-factor suppression) is kinematically accessible everywhere in the Sun and 
in fact enjoys an enhancement in its associated capture rate because of the reduced energy 
transfer and the related softening of form-factor suppression! This enhancement is by and 
large sufficient to compensate for the loss of helium (and partly oxygen) as a scatterer. 

4. Results 

In this section we present the numerical results for the capture rate of WIMPs in the Sun. 
Throughout we will assume a WIMP-nucleon cross-section of a n = 10~ 40 cm 2 , and note 
that the resulting capture rate scales linearly with the cross-section 8 . Using the DarkSUSY 
software package version 5.0 [27] we also give the associated muon yield on the Earth from 
different annihilation channels of WIMPs in the Sun. Many annihilation channels are in fact 
already constrained by present bounds on the neutrino flux from the Sun as we discuss below. 
Finally, we discuss the reach of future neutrino telescope experiments. 

In Fig. 3 we plot the capture rate in the Sun as a function of the WIMP's mass for 
two different choices of the most probable WIMP velocity in the halo. In the same figure 
we also depict the capture rate's dependence on the excitation energy, S. As 5 increases, 
the exponential suppression due to the nuclear form factor is curbed since less momentum 
is transfered. This results in a rapid rise of the capture rate as illustrated in Eq. A-14 in 
appendix A. However, when 5 becomes too large the capture rate diminishes rapidly because 
fewer shells in the Sun can participate in the capture. 

There are existing limits on the flux of muon-neutrinos from WIMP annihilation in 
the Sun from both underground detectors (BAKSAN [28], Super-Kamiokande [29], and 
MACRO [30]) as well as dedicated neutrino telescopes (AMANDA [31], BAIKAL [32]). The 
strongest bounds acually come from the Super-K results which we show in Fig. 4 alongside 
the muon yield for several annihilation channels plotted against the WIMP's mass. These 
channels are excluded by several orders of magnitude. On the right of Fig. 4, we plot the cor- 
responding limits on the branching ratios of the different annihilation channels. The bounds 
on direct annihilation into neutrinos deteriorate at higher WIMP mass because the more 
energetic neutrinos are further attenuated by matter in the Sun. 

Future neutrino telescopes (Antares [33] and IceCube [34]) are expected to have larger 
exposures and may provide even stronger constraints on WIMPs annihilation in the Sun. In 
Fig. 6 we show the expected reach for both hard and soft spectra together with the expected 
yield from the least constrained annihilation channels (cc and bb, both resulting in a soft 

8 cr n = 1CP 40 cm 2 is the typical WIMP-nucleon cross-section used in iDM models for fitting the DAMA 
results, see for example Ref. [3]. It is in fact somewhat of a lower bound since the needed cross-sections are 
almost always higher except when m x ~ 100 GeV and S < 50 keV [17]. 
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Figure 3: On the left pane we plot the capture rate against the WIMP mass. The solid curve 
correspond to an inelastic model with 5 = 125 keV and u© = 220 km/s whereas the dotted curve to 
vq = 254 km/s. In both cases we take a n = 10 -40 cm 2 . On the right pane we depict the growth of the 
capture rate as a function of the inelasticity for v© = 220km/s (solid) and u© = 254km/s (dashed). 
The upper two curves correspond to m x = 200 GeV (blue) and the lower two curves to m x = 400 GeV 
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Figure 4: On the left is a plot of the muon yield in the inelastic case (6 = 125 keV, a n = 10~ 40 cm 2 ) 
for different annihilation channels from top to bottom on the left: v^v^, t + t~ . Z a Z° , W + W~, ti, bb, 
and cc. The area above the thick (violet) curve is excluded by Super-K. On the right pane we plot 
the corresponding bound on the annihilation branching ratio for the respective channels against the 
WIMP's mass. 



spectrum). We assumed a muon threshold of 5 GeV which is appropriate for Antares, but 
probably too low for IceCube. 

Recently, motivated by the positron excess seen in PAMELA, it was suggested that DM 
may annihilate into two light bosons which are weakly mixed with the SM and so subsequently 
decay into light particles such as electrons, muons and pions [7]. As mentioned above, such 
channels cannot be probed by neutrino telescopes since those final states will stop in the Sun 
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Figure 5: Same as Fig. 4, but with Vq = 254 km/s. 



before decaying to neutrinos. However, gauge invariance requires these light bosons to also 
mix with the Z° and so have a small coupling to neutrinos as well. In the case of kinetic 
mixing the coupling to neutrinos compared to the coupling to charge is suppressed by m 2 /M|, 
where m is the mass of the dark vector-boson. Considering the enormous yield from direct 
neutrino production this might be observable when m > 3 GeV. However, in that case, the 
unsuppressed decay into charm quarks is likely to give a better bound. 




Figure 6: The bottom (top) solid curve shows the muon yield from WIMP annihilation in the Sun 
into cc (bb) which result in a soft muon spectrum for vq = 220 km/s ( vq = 254 km/s) on the left 
(right) pane. Antares' reach is shown for hard (soft) spectrum in the violet solide (dotted) curve. 
IceCube's reach is shown with the blue solid (dotted) curve for hard (soft) spectrum. 

5. Conclusions 

The iDM scenario offers the exciting possibility of extra structure in the DM sector which 
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explains the discrepancy between the DAMA results and limits coming from other direct 
detection experiments. It may also be responsible for some unexplained background events 
in those detectors employing heavier targets (XENON [35] and CRESST [36]). In this paper 
we showed that such a mechanism will result in a very large capture rate of WIMPs in the 
Sun. Current bounds on the neutrino flux from Super-Kamionkande strongly constrain the 
annihilation channels of such WIMPs into SM particles. Assuming a WIMP-nucleon cross- 
section of a n = 10 -40 cm 2 , it requires the branching ratio of annihilations into W + W~ , Z°Z°, 
t + t~, ti and neutrinos to be < 1%. Annihilations into bb and cc, which result in a softer 
spectrum, are less constrained. Future limits from Antares or IceCube may reduce these 
branching ratios by an additional order of magnitude. 

It is important to realize that these results are likely only robust to within ~ 50% due 
to the uncertainties in some of the input parameters (such as WIMP velocity distribution, 
iron distribution in the Sun, precise form-factors and etc.). Therefore, while the harder 
annihilation channels are strongly constrained, the softer ones are less so. Finally, it was 
recently shown in Ref. [37] that a dark disc (as compared to the spherical halo assumed in 
this work) would enhance the capture rate by up to an order of magnitude. If this effect is 
taken into account, even the softer annihilation channels are strongly constrained. 

Interestingly enough, the recent framework suggested in Ref. [7] unabashedly escapes 
these bounds. Motivated by the recent results from PAMELA, the DM in this picture an- 
nihilates predominantly into light leptons or pions through intermediate light bosons. Such 
light charged particles stop in the Sun before they decay and do not result in any observable 
neutrinos. Nevertheless, the bounds discussed in this paper restrict the coupling of DM to 
the SM. Aside from restricting the possible annihilation channels, it also has a direct impact 
on the possible mechanisms for the mediation of supersymmetry breaking to the extended 
DM sector. 

Note Added: We would like to thank the authors of Ref. [38] for the useful communi- 
cation regarding their results which are in qualitative agreement with ours. While they also 
consider the inelastic capture of WIMPs in the Sun, some of their input parameters are dif- 
ferent and serve to emphasize the inherent uncertainties in the computation when compared 
with our results. 

Acknowledgments: We would like to thank N. Arkani-Hamed, D. Lynden-Bell, A. 
Pierce, and N. Weiner for very useful discussions. I. Y. would like to thank the CCPP for 
their hospitality while this work was completed. L.-T. W. and I. Y. are supported by the 
NSF grant PHY-0756966 and the DOE grant DE-FG02-90ER40542. 

A. Capture Rate Computation 

In this appendix we summarize the steps, approximations, and figures used to compute the 
capture rate of inelastic DM in the Sun. We also provide an analytic approximation which 
help elucidate some of the features present in the numerical results. 
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We consider a WIMP with velocity u at infinity which is scattered in a region with escape 
velocity v. Hence, its total velocity in that shell is w 2 = u 2 + v 2 . In order to be captured, it 
must scatter down below v. The capture rate per unit shell in the Sun is given by [5], 

= J du M wn{w) (A-i) 

where f(u) is the WIMP's velocity distribution at infinity. fl(w) is the rate per unit time at 
which a WIMP of velocity w scatter to a velocity less than v. It is given by, 

n/ \ / \ Qrnax ~ Qcap / a n\ 

n(w) = (naw) — (A-2) 

where the first factor, naw is just the rate of scattering, and the second factor embodies the 
probability of capture. We now address each component separately. 
The inelastic non-relatvisitic cross-section is given by, 

f 2 Z 2 + f n (A-Z) 2 \ 

JjL ji -J a n (A-3) 

with a ne being the WIMP-nucleon reduced mass, and a n is the WIMP-nucleon elastic cross- 
section. fp( n ) determines the relative contribution of protons (neutrons) and we normalize 
our results to /„ = f n = 1. We choose a fiducial value of a n = 10 -40 cm 2 . The final rate is 
easily scalable with these quantities. 

The second component in Eq. (A-2) contains Q m ax {Qmin) which is the maximum (min- 
imum) energy transfer possible in the scattering process. Q cap is the minimal energy transfer 
needed for capture, 



Qmax = ^m x w 2 1-Ar 1-— Jl ^: \-6 (A-4) 





Q imn = ^m x w z 1-^V 1 + — Wl \-S (A-5) 



Qcap = T, m x { w2 -V 2 ) - 5 (A-6) 



where u is the WIMP-nucieus reduced mass, and 5 is the inelastic energy gap. Including 
nuclear form factors, the formula for the capture rate per unit time gets modified to, 

Wrnax Wmin JQ cap 

where, 

«°=2^ fl = 10 " 3 ( M1 (w) +M ) ™ <A - 8) 
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We assume that the WIMP's velocity distribution in the halo is given by a Maxwell- 
Boltzmann distribution as seen by an observer moving with velocity v & 

f(x)dx = ^ * ^^^ sinh^,) U x (A _ Q) 

Tfl x -y 7T ZXTj 



where p = 0.3 GeV cm 3 is the local mass density and the dimensionless quantities x and r\ 
are given by, 



u 2 



x 2 = (A-10) 



= $ (A-ll) 

with vq = 220km/s. 

It is possible to obtain approximate analytic results for the capture rate as follows. Be- 
cause of the phase-space suppression in Eq. (A-3) the scattering has larger support when 
5 < pw 2 /2. Hence we can Taylor expand the radical in the expression for the momentum 
transfer and get, 

Q max = \™A^ (A 2 + (3_x 2 - (^hi+^A *_\ _ S (A-12) 
2 /?_ \ \m x - m N J m x v z /2 ) 

with f3± = 4m x m N /(m x ±m N ) 2 , A 2 = v 2 (3-/vq (our notation here follows that of Gould [5]). 
In this case, ignoring form-factor for the moment, the integral over the velocity distribu- 
tion can be done exactly and the results are identical to those found in Ref. [5] with the 
identification, 



A 2 ^A 2 -[ •^LZ^IL (a-13) 



m x + m N 

m x — m N J m x v 2 /2 

This tends to reduce the capture since it diminishes the maximum energy transfer and hence 
the probability for capture. Notice that the additive term involving 5 in the energy transfer 
drops out since in the absence of form-factors the capture rate involves only differences in 
energies. 

This term becomes significant when including the form-factor as it serves to increase the 
capture rate. Using Eq. (A-13) it is possible to obtain an approximate analytical expression 
in this case as well. Regardless, the most important feature is that the capture rate grows 
exponentially with S, 

dC _ 2anp x v s/Qo 

dV ^m x 1 ' 

where the ellipsis denote the usual expression involving the form-factor corrected by 5. This 
behavior reflects the fact that as the inelastic threshold increases, less energy is transferred 
to the nucleus and the form-factor suppression is correspondingly smaller! Therefore, in the 
case of inelastic transitions, the reduction due to form factor is not as severe. In particular it 
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allows iron to assume the role of the prime contributer to capture, a role it would have played 
in the elastic case as well if not for the form-factor suppression. 

All together, inelasticity enters in two significant ways. First, it tends to reduce the 
capture rate because of the diminishing phase-space. Second, since less energy is transferred 
to the nucleus, the reaction is more coherent, form- factor effects are diminished, and the 
capture rate enjoys an increase as compared with the elastic case. For iron in particular, the 
second effect is much more important as discussed in the text above. 

B. Density Computation 

In order to compute the density of the WIMPs accumulated in the Sun we approximated 
the Sun's potential with an analytic formula which allows for an exact solution for the orbits 
[22, 23]. It interpolates between an harmonic potential close to the core and an inverse fall-off 
at long distances, 



(in here the potential energy is computed with respect to the origin where it is zero rather 
than at infinity). A fit to the Sun's actual potential [25], gives b = 0.0884-Rq and yields an 
approximation which is roughly 10% — 20% accurate. The orbits are found by first expressing 
the potential U in terms of r, which by an adequate choice of units can be written simply as, 



The dependence of the azimuth and the time variable on the potential (and hence on r) can be 
found by usage of the constants of motion and a straightforward integration. The interested 
reader can find the details of the solution in Ref. [23] . 




Figure 7: A typical event in which a WIMP approaches the Sun, collides with an Iron nuclei and 
gets captured in a bounded rosette- like orbit. The solid-blue line is the approach path from infinity 
to the point of collision. The dashed-green line is the orbit after collision which is normally contained 
within the Sun that is delineated by the circle. 

To find the final density, we incorporated the exact solution for the orbits into a Monte- 
Carlo simulation. We let the WIMP's start at infinity with an impact parameter which is 
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uniformly distributed and a shifted Maxwel-Boltzmann velocity distribution with the rota- 
tional speed of the local standard of rest vq = 220 km/s and v & = vq. The WIMP is then 
allowed to scatter at any radius in the Sun, assuming its orbit reaches that point and its kinetic 
energy is sufficiently large to overcome the inelasticity. We take the probability of scattering 
as a function of radius to be proportional to the local density divided by the WIMP's velocity 
at that radius. The scattering is assumed to be isotropic in the center of mass frame, and the 
kinetic energy of the outgoing WIMP is restricted to be below the capture rate. After the 
collision, the kinetic energy of the outgoing WIMP and its angle with respect to its original 
direction of motion are used to determine its new orbit. Fig. 7 depicts a typical approach 
and capture of a WIMP by the Sun. If any further collisions are energetically possible then 
they are executed following similar steps as described above. Once no further collisions are 
possible the orbit is recorded and stored. 

The density is computed by simulating 10,000 capture events and then randomly sampling 
the radii of the resulting orbits. It is important to sample uniformly in time rather than radius 
because the particle actually spends most of its time close to the apogee. 
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